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Mesons at finite baryon density in (2 + l)d 

Costas G. Strouthos 
Department of Physics, Duke University, Durham, NC 27708-0305, USA 

We discuss the critical properies of the three-dimensional Gross-Neveu model at nonzero 

■1 temperature and nonzero chemical potential. We also present numerical and analytical 

, results for the in-medium interaction due to scalar meson exchange. Further, we discuss 

' in-medium modifications of mesonic dispersion relations and wavefunctions. 
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Chiral phase transitions and the spectrum of excitations in strongly interacting 
, matter at finite baryon chemical potential remain interesting challenges. Strongly in- 

teracting systems are intrinsically non-perturbative and therefore most of our knowl- 

■ edge about the relevant phenomena comes from lattice simulations. Unfortunately, 
, the complex nature of the determinant of the QCD Dirac operator at finite chemical 

potential makes it impossible to use standard simulation techniques to study Fermi 

■ surface phenomena in Euclidean simulations. 

OO , In order to understand what ingredients might play a decisive role in more 

^ I complex systems such as gauge theories, we have studied the simplest non-trivial 

' model simulable with ^ using standard algorithms, namely the three-dimensional 

Gross-Neveu model (GNM3). Its Lagrangian in Euclidean metric is written in terms 
of 4A'^y-component spinors ■0, ip as 



Oh! 2 

C = i;{i^ +m)^l^-^i^l,^Pf. (M) 
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In the chiral limit m = the model has a global Z2 symmetry ip 1— > 75'0, ip —^75- 
;J] ■ At tree level, the field a has no dynamics; it is trully an auxiliary field. However, 

^ it acquires dynamical content by dint of quantum effects arising from integrating 

out the fermions. The model is renormalizable in the l/Nf expansion unlike in the 
loop expansion.^ Apart from the obvious numerical advantages of working with a 
relatively simple model there are several other motivations for studying this model. 
At T = ^ = for sufficiently strong coupling g^, chiral symmetry is spontaneously 
broken by a condensate (tpip) 7^ leading to a dynamically generated fermion mass 
gap given by Mf = g'^{'ip'tp) ^ m in the large- A'^y approximation. The spectrum of 
excitations contains both baryons and mesons, i.e. the elementary fermions / and 
the composite // states. The critical coupling g'^ at which the gap Mf/Auv 0, 
defines an ultra-violet stable fixed point of the renormalisation group at which an 
interacting continuum limit may be taken. This picture has been verified at next- 
to-leading order in 1/A^/ and by Monte Carlo simulations with finite 

In the next section we will discuss the issue of chiral symmetry restoration at 
nonzero fermion chemical potential. In section |31 we will discuss Fermi surface phe- 
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nomena, such as long range interparticle interaction in medium. We will also present 
results for dispersion relations (section ^ and wavefunctions (section ISJ of massless 
particle-hole excitations on or near the Fermi surface. A study of the dispersion 
relation of spin-^ excitations around the Fermi surface has shown that there is no 
mass gap characteristic of a BCS instabilit}® and the values of the Fermi liquid pa- 
rameters (Fermi momentum, Fermi velocity) are in good agreement with the 1/-/V/ 
predictions. However, in a different study of a related model, the {2 + \)d NJL model, 
which has an SU{2) (8) SU{2) chiral symmetry, it was shown that although the BCS 
instability is also absent the properies of the quasi-particle spectrum is determined 
by physics outside the scope of the 1/Nf expansion.® 

§2. Chiral phase transitions 

The action of the NJL model remains real even after the introduction of nonzero 
chemical potential fi, which means we can study the physics of the high density 
regime using standard Monte Carlo techniques. In the presence of a Fermi surface 
with Fermi momentum pp the creation of // pairs with zero net momentum is 
energetically suppressed, because one can only create particles with p > pp- So as 
the fermion number density ry(/x) grows the chiral symmetry breaking is suppressed. 
The large- A'^^ description of the jU 7^ chiral phase transition predicts a first order 
transition for T = and a continuous transition for T > 0^ Furthermore, the 
critical value of the chemical potential fic is equal to the value of the fermion mass at 
/i = 0, which indicates that materialization of the fermion itself drives the symmetry 
restoration transition. Interactions as expected decrease Hc below the mean field 
result.® Work by Stephano'^^^ suggests that any non-zero density simulation which 
incorporates a real path integral measure proportional to det (MM^^) is doomed to 
failure due to the formation of a light baryonic pion from a quark q and a conjugate 
quark q'^. The NJL model however, does not exhibit such a pathology, because 
the realization of the Goldstone mechanism in this model is fundamentally different 
from that in QCD. In the NJL model the Goldstone mechanism is realized by a 
pseudoscalar channel pole formed from disconnected diagrams and the connected 
diagrams yield a bound state of mass ~ 2m j.^' This implies the absence of a light 
qq'^ state. 

As expected our simulations of GNM3 with Nj = ^ did not provide any evi- 
dence for the existence of a nuclear liquid- gas transition at // < /ic- It was shown 
jjjioj i^j^at in the standard four-fermi models there is no saturation density for stable 
matter. In order to get the saturation features the authors oPJ introduced a chirally 
invariant scalar-vector interaction term which cures the binding problem. Further- 
more, our results showed that the second order nature of the T 7^ 0, = transition 
remains stable down to low T and large /x. In Fig^ (left) we plot the normalized 
order parameter E(T, fi)/ U{T,0) {U = (a)) as a function of fj,/mf at different val- 
ues of T. It is clear from the shapes of these curves that the transition becomes 
sharper as we decrease the temperature. In Fig^ (right) we plot the normalized 
fermion number density as a function of the chemical potential at different values 
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Fig. 1. Chiral order parameter (left) and fermion number density (right) vs. chemical potential at 
various values of T. 



of T. In the limit T ^ we see that the fermion density is strongly suppressed 
before the transition and then jumps discontinuously. By performing detailed finite 
size scaling analysis which allowed us to distinguish between second order and weak 
first order transitions we showed that the tricritical point lies on the section of the 
phase boundary defined by T/T^ < 0.23, > 0.97.1^ It was also shown some time 
agJS that the second order transitions belong to the 2d Ising universality class in 
accordance with the dimensional reduction scenario. 



§3. The sigma propagator in medium 



In this section we examine the auxiliary scalar propagator D^j in the presence 
of non-zero baryon density. Although at tree level a is non-propagating, the leading 
order l/Nf expansion at = predicts that it acquires dynamics through quantum 
corrections due to virtual qq pairs, resulting in a propagator of the forirP''^ 

D,{k^) = 4 ^ ' , ^ :r-. T^- (3-1) 

5^ 2r(2 - f ) + 4M|)F(1, 2 - f ; I; 

Immediately we see the difference between this model and QCD. For k'^ <^ Mj 
the hypergeometric function F ~ 1, implying that to this order the a resembles a 
weakly-bound meson of mass Mu = 2Mj; however, the hypergeometric function in 
the denominator indicates a strongly interacting qq continuum immediately above 
the threshold 2Mj. This implies that if truly bound, its binding energy is 0{1/Nf) at 
best (to our knowledge there have so far been no analytic calculations), implying little 
if any separation between pole and threshold. A recent study using the maximum 
entropy method has shown evidence for a nonzero binding energy in the a channel for 
finite Nf = 4.^3 Since all residual interactions are subleading in l/Nj, we surmise 
that all other mesons are similarly weakly bound states of massive fermions, and 
hence effectively described by a two-dimensional "non-relativistic quark model" . 

Similarly, at nonzero chemical potential fi to leading order in the l/Nf expansion, 
we have 



D-\k) = l-n{k;fi) 
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(3-2) 



where 77 is the virtual fermion-antifermion vacuum polarization bubble. We have 
used the gap equation at /i = to express in terms of the zero density gap 
1^0 and assumed that the gap vanished in the integral defining 77 for fi > Hc- The 
functional form of the propagator in the Hard Dense Loop (HDL) approximation 
(/cq, 1^1 *C /u) in d + 1 dimensions is given 



D-Hko,k) 



g2^d 3 



(4vr)2r(f)(3-d) 



1 



(3-rf)(fcg + |fc| 



ko\k\^ 



(3-3) 
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4^ 



1 + 



[ko + {ki + m-^] 



An interesting observation is that in the static limit fco = the momentum depen- 
dence of eq. (|3-4j) vanished for the value d = 2 implying complete screening for the 
static potential due to a exchange. Furthermore, for /c = 0^ is proportional to a 
conventional zero-momentum boson propagator with mass 



2\/Ai(/^-^c), 



(3-4) 



implying that just above the transition there is a tightly bound state. In Fig. I^Jwe 
show the results for the sigma correlator obtained from numerical simulations with 
Nf = 4 and l/g^ = 0.75 on 48 x 32^ lattices. The value of the critical chemical 
potential = 0.17(1). The numerical results for M^- fall 20 — 30% lower than the 




Fig. 2. The zero momentum o timeslice correlator for 4 values of /i (left) and the Ma as a function 
of ji (right); the line is the leading order 1/Nf prediction. 

leading order prediction, which is consistent with a correction of 0{Ny^). 

§4. Mesonic dispersion relations 

In this section we investigate mesonic states by measuring the connected con- 
tribution to the correlation functions Cr{k,xo) = X^^(ir(0, 0)j|,(x, xo))e~*'^'^ where 
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the bilinears jr{x) are defined with scalar, pseudoscalar or vector quantum numbers. 
In terms of staggered fermion fields X the operators are 

h{x) = XxXx] j-ysi^) = ^xXxXx] h{x) = ^[XxXx+i + Xx+iXx], (4-1) 

where rjix = (—1)^", t/2x = (— g^j^j _ (•_]^-)a;o+xi+2:2^ bgfQj-g -^g study 
the timeslice correlators in each channel as a function of spatial momentum /c || x. 
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Fig. 3. Pseudoscalar correlator C75(|fc|,a;o) at 4 different /i for momenta \k\ — (filled circles), 
(filled squares), j (filled diamonds), |^ (filled up triangles), ^ (filled down triangles), |^ (open 
circles), ^ (open squares), ^ (open diamonds) and ^ (open up triangles). 

It was showiJ^ that in free field theory for fi > fic the the xq — > limit of Cr(A:, xq) is 
dominated by a continuum of particle-hole pairs at or near the Fermi surface, which 
effectively cost zero energy to excite. The generic result is that for \k\ < 2/x the 
decay is algebraic, with 



Cr{k,xo) oc 



\k\ < fi 
\k\ ~ 2//. 



(4-2) 



Only once \k\ > 2jjL does it become kinematically impossible to excite a pair with 
zero energy, resulting in exponential decay: 



Cr{k, Xq) oc Xq ^ exp( — — 2//)xo j . 



(4-3) 



The sequence of plots in Fig. I^Jis in qualitative agreement with these findings. For 
each /i there is a particular value of \k\, highlighted with a solid line in the plots, for 
which the temporal falloff is particularly slow, corresponding to \k\ ~ 2^ (e.g. for 
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= 0.4 the slow falloff occurs for |fc| = ^ ~ 0.785). For \k\ larger than this value 
the decay is much steeper, although only for = 0.2 does it resemble an exponential 
form. Because of the technical difficulties in treating correlators with power-law 
decays in a finite volume we have made no attempt to fit the numerical data for 
Cy^{\k\,xo) to a functional form. 

Light states in the p channel have long been of interest. Here we discuss tentative 
results in the vector channel, corresponding to the quantum numbers of the p meson. 
In this case for \k\ > we can distinguish between C-^n (k), in which the component of 

the vector is parallel to k, and C^j_{k). In Fig. |l]we plot the correlator for several \k\ 
values at ^ = 0.6 in each case. For C^^ the curves are qualitatively very similar to 
those of Fig. 131 at p = 0.6, with a distinguished momentum \k\ = The correlator 
C-^ii is much smaller in magnitude, and is consistent with exponential rather than 

algebraic decay. The lines are fits of the form C-y^^ {\k\, xq) = A{e~^^° + g-^C^t-^o)). 

For small \k\ the resulting resembles that of a massless pole, ie. E = /3o|^|; 

with velocity (5q ~ 0.5. It should be mentioned, however, that although and 
still differed, no evidence for a massless pole was seen in the data at p = 0.8. 
A potential problem for this picture is that the extracted velocity /?o ~ 0.5 is less 
than /^iT' ~ 1, implying that quasiparticles would experience damping via Cerenkov 
radiation of zero sound. Clearly further work exploring the systematic effects of 
varying p, g^, k and volume will be needed for a more complete understanding to 
emerge. 
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§5. Mesonic wavefunctions and Friedel oscillations 

In this section we expose another characteristic of the Fermi surface by study- 
ing the spatial correlations between fermions and antifermions in various mesonic 
channels, probed via the wavefunction iP'(x) defined by 

^r{x) = lim a^^i(x = 0) V(gg(0,0;y,xo)rg5(0,0;y + x,2;o)r), (5-1) 

y 

where as usual F projects out the quantum numbers of the channel of interest. In 
the GNM3 and related models ^{x) is technically much easier to define and measure 
than in QCD-like theories where they are gauge-dependent. Meson wavefunctions in 
the GNM3 model have been studied at T, ^ = in,-^' where further technical details 
are given. 

It is easy to show^- that in the large- A'^y approximation the PS wavefunction in 
the chirally broken phase is given by 

hm C(y;xo) ~ -e'^^- exp • (5-2) 

The general profile of the wavefunction is a gaussian with width increasing as ^JxQ. 
In the chirally symmetric phase we get, 

hm C(y;xo) ~ ^e-^'^-" Jo(/z|y|). (5-3) 

a;o->0 Xq 

The wavefunction profile no longer changes with xq, but instead oscillates with a 
spatial frequency determined by /i, which to this order may be identified with the 
Fermi momentum kp. The oscillations observed in C{y;xo) are characteristic of a 
sharp Fermi surface and are reminiscent of oscillations in either the density-density 
correlation function, or the screened inter-particle potential, in degenerate systems 
generically known as Friedel Oscillations^^ In Fig. I^we plot the PS wavefunctions 
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Fig. 5. Scalar wavefunction at jj. — 0.4 (left) and /j, — 0.6 (right). 



at = 0.75 for ^ = 0.4,0.6. The simulation data are connected by solid lines. 
As fi, and hence kp increases the oscillations decrease in wavelength in accordance 
with the theoretical expectation and provide a graphic illustration of the presence 
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of a sharp Fermi surface. The dashed hnes show measurements taken with the 
same lattice parameters but with the interaction switched off. The disparity with 
the interacting theory is small, though increasing with showing that the free 
field description of the oscillations is qualitatively correct. A possible explanation 
is the infinite value predicted for the Debye mass in section |21 implying that inter- 
actions between static quarks at non-zero \x\ are completely screened. Also shown 
is the theoretical form (eq. I5v-i|) ^[x) = jQ{kpx) with kp = ^, showing good agree- 
ment with the data for small x ^ 2kp^ . Unfortunately it appears hard to obtain 
more quantitative information, such as an independent fit for the Fermi momen- 
tum kp, because \Jo{kx)\ decays only as x~2. This means that fits should not only 
include the backwards-propagating signal Jo{kp{Ls — x)) but also image contribu- 
tions J{){kp{nLs — x)) - our attempts to find a satisfactory fit were unsuccessful. 
The figures therefore simply show both the "forwards" Jo(/ix) and "forwards-and- 
backwards" Jo(/Ux) -|- Jo(/i(Ls — x)) forms, showing that neither gives a satisfactory 
description of the data over the full range of x and /i. We also checked that there 
are no significant differences among the various mesonic channels, implying that in 
contrast to the situation at = 0,^ effects due to eg. a exchange are very hard to 
detect. 

§6. Summary 

The four-fermion models remain the only interacting fermionic field theories 
both simulable by standard lattice methods at /x 7^ and displaying a Fermi surface, 
thought to be an essential figure of dense quark matter. We have shown that the 
/X 7^ 0, T = transition is strongly first order and a tricritical point exists in the 
(/i, T) plane at very low T. 

We also discussed in-medium effects on the character of mesons. The large- A^j 
auxiliary boson propagator has a branch-cut in the complex-A; plane which is modified 
to an isolated pole in the high density phase. Lattice simulations with Nf = 4 have 
verified this with reasonable precision. We have also shown that the connected (flavor 
non-singlet) diagrams instead of showing exponential falloff with Euclidean time, 
they generically decay algebraically, signalling the presence of massless particle-hole 
excitations. We also presented tentative evidence for a massless pole in the vector 
channel, which is possibly a manifestation of zero sound. Furthermore, a graphic 
confirmation of the presence of a sharp Fermi surface was shown by the oscillatory 
behavior in mesonic wavefunctions, which resemble the Friedel oscillations in many- 
body physics. 
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